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Abstract 



The modeling of UV and optical spectra emitted from the symbiotic system AG Draconis, adopting collision of the winds, predicts 



(N 



soft X-ray bremsstrahlung from nebulae downstream of the reverse shock with velocities > 150 kms and intensities comparable 



Qh to those of the white dwarf black body flux. At outbursts, the envelop of debris, which corresponds to the nebula downstream of the 
^ high velocity shocks (700- 1 000 km s~ 1 ) accompanying the blast wave, absorbs the black body soft X-ray flux from the white dwarf, 
explains the broad component of the H and He lines, and leads to low optical-UV-X-ray continuum fluxes. The high optical-UV 
flux observed at the outbursts is explained by bremsstrahlung downstream of the reverse shock between the stars. The depletion of 
1—1 C, N, O, and Mg relative to H indicates that they are trapped into dust grains and/or into diatomic molecules, suggesting that the 

, , collision of the wind from the white dwarf with the dusty shells, ejected from the red giant with a ~1 year periodicity, leads to the 

fV| U-band fluctuations during the major bursts. 
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1. Introduction 

Circumstellar and circumbinary nebulae throughout symbi- 
otic systems (SS) emit line spectra and continuum in the ra- 
dio, infrared (IR), optical, UV, and X-ray range (Angeloni et al. 
2010 and references therein). 

The two main SS components, generally a white dwarf (WD) 
and a red giant (RG) rotating in bound orbits, lead to character- 
istic time periodicities. Moreover, they dominate the continuum 
spectral energy distribution (SED) in the UV-soft X-rays and 
IR-optical ranges, respectively. 

The RG flux is well reproduced by a black body (bb). The 
analysis of symbiotic systems shows that the flux from the hot 
component star is approximated by a bb which peaks in the 
far UV. It was rarely observed, beyond a few data in the bb 
low frequency tail. Consequently, the temperature of the hot 
star was obtained indirectly from the characteristic spectral line 
ratios emitted from the nebulae (Angeloni et al. 2010). 

AG Draconis (AG Dra) is classified as a yellow symbiotic 
binary containing a K2III giant (Miirset & Schmid 1999) and a 
WD accreting from the giant's wind on a 549 day orbit (Fekel et 
al. 2000). The orbital inclination is relatively low with i ~30°- 
45° (Mikolajewska et al. 1995) or i ~ 60+ 8.2° (Schmid & 
Schild 1997). In fact, during the orbital motion the WD is not 
eclipsed. 

The light curve of AG Dra shows numerous bursts with am- 
plitude spanning 1-3 mag in U. Gonzalez-Riestra et al. (1999, 
hereafter GVIG99) identified "cool" and "hot" outbursts differ- 
ing in Zanstra temperature and light curve profile (Skopal et al. 
2009). Tomov & Tomova (2002 and references therein) report 
from the analysis of the photometric data that the visual light 



growth during the optical brightening is due to increased radia- 
tion of the circumbinary nebula. Thermonuclear events on the 
surface of the WD (Mikolajewska et al al. 1995) and expansion 
and cooling of the compact secondary (Greiner et al. 1997) 
were also conjectured to explain the outbursts. 

Besides the characteristic outbursts in the U band, AG Dra 
is the brightest system at soft X-ray energies of all SS. Viotti 
et al. (1998) monitored the super soft X-rays from the AG Dra 
source by EXOS AT and ROS AT and the UV flux by IUE. It was 
found that during all the monitored outbursts the visual magni- 
tude was associated with an increase in the UV flux, while the 
X-ray count rate largely faded, namely the count rate dramati- 
cally dropped at the time of the 1985 and 1986 maxima (Viotti 
et al. 1998) and between 1992 and 1996 (Greiner et al. 1997). 
Actually, Viotti et al. (1998) discovered the systematic opti- 
cal/UV and X-ray flux anti-correlation, which was confirmed 
by Skopal et al. (2009), Gonzalez-Riestra et al. (2008), etc. 
observations. 

Analyzing new data by XMM Newton and FUSE, Skopal 
et al. (2009) recently concluded that the soft X-ray emission 
is produced by the WD photosphere. They claim that X-ray 
and far UV observations determine the bb temperature unam- 
biguously. Modeling the UV and super-soft X-ray data in AG 
Dra, they concluded that the observed anti-correlation is caused 
by the variable wind from the hot star, explaining that the en- 
hanced hot star wind gives rise to the optical bursts by repro- 
cessing high energy photons from the Lyman continuum to the 
optical/UV. 

This paper is motivated by the analysis of AG Dra XMM- 
Newton observations (Jansen et al. 2001) performed during the 
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2003 hot burst and during quiescent phases by Skopal et al. 
(2009), who found that the soft X-ray spectrum at quiescence 
corresponds to bb radiation. We would like to investigate by 
consistent modeling of the line and continuum spectra whether 
the bremsstrahlung downstream of shocked nebulae might con- 
tribute to the soft X-ray emission at different phases. 

Significant information about the physical conditions of the 
SS components can be obtained by line spectrum analysis. The 
intensities and profiles of some spectral lines emitted from AG 
Dra at different epochs were analyzed in detail (e.g. the Ha and 
Hell 4713 lines by Tomov & Tomova 2002, the OVI 1032,1035 
doublet and Hell 1640 line by Skopal et al. 2009). The most 
comprehensive and complete sample of spectra in number of 
lines at different close phases in the optical and UV ranges was 
provided by GVIG99. Both the line and continuum spectra 
were explained by photoionization of circumstellar matter, e.g. 
the wind from the RG, and of circumbinary matter by the hot 
star radiation flux. 

Viotti et al. (1989) noticed, however, that in some models 
X-ray are interpreted as bremsstrahlung emission from the col- 
lision of the winds from both components (e.g. Kwoks & Leahy 
1984, Girard & Willson 1987). The steady and intense X-ray 
emission of AG Dra during quiescence indicates the presence 
of an effective heating process. 

Therefore in this paper we calculate the spectra using com- 
posite models which account consistently for shocks and pho- 
toionization. We will follow the method adopted by Angeloni 
et al. (2010) for other symbiotic systems, cross-checking the 
calculations of the continuum by those of the line ratios. We 
will focus on the optical-UV-X-ray domain, leaving the radio- 
IR range to the next investigation of the AG Dra SS. 

The models presented by Skopal et al. (2009) for AG Dra are 
valid on a large scale, namely the continuum emitted from the 
nebulae was not constrained by the analysis of the line spectra. 
The collision of the winds cannot be neglected because it leads 
to high temperatures of the gas. In this paper we would like 
therefore to investigate whether, adopting models accounting 
for both photoionization and shocks, it is possible to add some 
information to the interpretation of AG Dra spectral features, in 
particular, and of SSs, in general. 

The SS modeling and the calculation of the spectra are ex- 
plained in Sect. 2. The spectra are analyzed in Sect. 3. Results 
are presented in Sect. 4. Discussion and concluding remarks 
follow in Sect. 5. 



2. The models 

Alternating episodes of accretion and wind collision were ev- 
ident during the evolution of symbiotic systems, e.g. CH Cyg. 
The epochs treated in this work on AG Dra are close to out- 
bursts. We consider that the accretion disk was blown up and 
had no time to form. In fact, the observed spectra do not show 
the typical double peaked lines. So we adopt the colliding wind 
model of Girard & Willson (1987), which applies close to the 
orbital plane of the SS. 



2.1. Wind collision theory 

Girard & Willson (1987) explained that in a widely separated 
binary consisting of a late giant and a hot WD, the late star has 
a stellar wind with a mass loss rate of 10 7 -10 5 M yr~' and a 
velocity ~10-30kms~'. The hot star will accrete hydrogen rich 
material from the late red giant wind. Theoretical models by 
Paczynski & Zytkow (1978) show that an accretion of hydrogen 
rich material at a rate of 10~ n to 10~ 10 M yr~' will generate 
periodic hydrogen shell flashes, with each outburst lasting on 
the order of decades. The increased flux from such an eruption, 
particularly in the UV region, and the accompanying increase 
in radiation pressure will tend to eject the hot star's hydrogen 
atmosphere in the form of a stellar wind with velocities of 1000- 
3000 km s" 1 (Conti 1978), which may last for decades. Both the 
winds from the WD and the red giant are supersonic, so shocks 
will form. 

Since the end of the 1990, it was suggested by the analysis of 
SS line spectra (Nussbaumer 2000 and references therein) that 
both stars lose winds which collide in the interbinary region and 
outward, creating a network of shock-fronts. A picture of the 
colliding wind model adapted to the presence of dust shells is 
given by Contini et al. (2009c, fig. 2). 

The scenario shows the collision of the winds from both stars 
leading to two shock-fronts which propagate between the stars: 
the strong one, dominating the spectrum, propagates in reverse 
towards the WD, and the weaker one, not always important, 
propagates towards the RG. Moreover, a shock-front expands 
out of the system through the circumbinary medium. 

2.2. Calculation of the spectra 

Composite models considering both photoionization and 
shocks in a consistent way were adopted (Contini 1997 - An- 
geloni et al. 2010). 

In the past years we carried out the analysis of SSs (WD + 
RG + nebulae + dusty shells) of both line and continuum spec- 
tra emitted from the nebulae downstream of shock-fronts and 
of dusty shells ejected by the RG. In particular, the observed 
SED was disentangled into the different components by adjust- 
ing model calculations to the data in the different frequency 
ranges. The models were constrained by cross checking the line 
ratios with the continuum. Comparison with the absolute ob- 
servational data determines the radius of the nebulae throughout 
the SS. 

The collision of the outburst high velocity ejecta with the 
RG wind and/or with the ISM surrounding the WD on the side 
opposite the RG (Contini et al. 2009a), induces the high tem- 
peratures downstream of the shock-fronts and, consequently, to 
X-ray emission (Angeloni et al. 2010, fig. 4), while in the UV 
domain the WD bb flux dominates. However, the data in the X- 
ray range were not always sufficient to distinguish the origin of 
the far-UV - soft X-ray emission. Both the bb tail from the WD 
and bremsstrahlung emitted from nebulae downstream of rel- 
atively fast shocks were proposed (e.g. Contini & Formiggini 
1999, fig.4). 

In particular, the nebulae downstream are photoionized and 
heated by the flux from the WD and collisionally by the shock. 
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The gas entering the shock-front is compressed and heated in 
the post shock-region. The flux from the WD and the secondary 
diffuse radiation flux are calculated throughout a large number 
of slabs downstream (up to 300) by radiation transfer. The gas 
cools down by free-free, free-bound and line emission. The 
temperature drop is calculated smoothly downstream through- 
out the slabs whose geometrical thickness is automatically cal- 
culated by the temperature gradient. The computation halts 
when the the gas is cool enough in the radiation bound case, 
or, in the matter bound case, when the edge opposite to the 
shock-front is reached. 

We use for the calculation of the spectra the code SUMaQ 
which simulates the physical conditions of an emitting gaseous 
nebula under the joint effect of photoionization from an external 
source and shocks, and in which line and continuum emission 
from gas are calculated consistently with dust reprocessed radi- 
ation as well as with grain heating and sputtering. 

3. The spectra 

In AG Dra, the broad and narrow components of the line 
profiles reveal high (700 to 3000 kms 1 ) and low (100-160 
km s _1 ) velocities at the epochs corresponding to the outbursts. 
The contemporary presence of high and low velocity shocks 
between the stars with a velocity gap by a factor > 10 in- 
dicates that the matter is highly fragmented by instabilities 
at the fluid interface (Richtmyer-Meschkov, Kelvin-Helmoltz, 
Rayleigh-Taylor) created by the outburst. Alternatively, the 
high velocities may characterize the shocks accompanying the 
blast wave created by the outburst in directions opposite to the 
RG. These shock-fronts can be treated by the dynamical evo- 
lution described by Chevalier (1982) in supernova atmospheres 
(Contini et al. 2009a). We will see in the following that the lat- 
ter hypothesis gives a more acceptable explanation of AG Dra 
observations. 

Following our method, first we model the line ratios, which 
constrain the model, then we cross check the results by com- 
paring the calculated continuum SED with the data. 

To constrain the models, the spectra should contain enough 
significant lines, such as those from the same element but from 
different ionization levels which give information about the 
conditions of the emitting gas independently of relative abun- 
dances, and lines from the same level of different elements 
which lead to different stratifications of the ions downstream, 
depending on the ionization potential. Comparison of line ra- 
tios in the UV and in the optical range also helps understanding 
the gas conditions within the SS. 

The SED reveals the relative importance of the different con- 
tributions to the continuum (the radiation from the stars, the 
bremsstrahlung from the nebulae, and the reprocessed radiation 
from dust grains). Generally, in SSs, the SED shows that dif- 
ferent nebulae contribute to the spectra. Their relative contribu- 
tions must be consistently accounted for also in the calculations 
of the line ratios (Contini et al. 2009b) leading to an iterative 
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calculation process which stops when the calculated spectra and 
the data are well tuned. 

The data are observed at Earth, while the calculations refer to 
the nebula. By adjusting the calculated to the observed absolute 
continuum by the factor rj~ r 2 Id 2 , which depends on the radius 
r of the nebula and on the distance d to Earth, some information 
about the location of the nebulae is obtained. 

3.1. UV and optical lines 

We refer to the UV and optical spectra observed by GVIG99 
between June 1979 and February 1996. To get a first estimate 
of the conditions in the nebulae at different epochs, we com- 
bined the spectra in the UV and in the optical range observed 
at adjacent days. This allows to refer also the UV line ratios 
to Hfi which is one of the most significant lines. The analy- 
sis of the line ratios reveals both the physical conditions and 
the relative abundances of the elements throughout the nebulae. 
We thus selected the UV spectra at 8 epochs: 4-11-1985 (JD 
2446374), 18-1-1986 (JD ..6449), 30-4-1986 (JD ..6551), 27-2- 
1993 (JD ..9046), 4-7-1994 (JD ..9538), 6-12-1994 (JD ..9693), 
28-7-1995 (JD ..9927), and 14-2-1996 (JD 245128) which are 
randomly chosen phases. The selected dates are indicated upon 
the diagram presented in Fig. 1 (bottom), where we report the 
profile of the Hj6 flux observed during the 1985-1986 and 1994- 
1996 active phases by GVIG99 (fig. 4, bottom diagram) and by 
Mikolajewska(1995). 

Optical line observations are given by GVIG99 (table 2). The 
dates closest to the 8 epochs chosen for the UV observations 
are: 25-10-1985,21-1-1986,22-5-1986, 10-1-1993, 1-7-1994, 
26-12-1994, 19-8-1995, and 3-1-1996, respectively. 

The observations (GVIG99) show that the spectral region in- 
cludes strong emission lines of hydrogen, helium, and weaker 
lines e.g. Fell and NIII. During the active phases the strongest 
emission lines, Hell 4686, H/3 and Hy, display extended wings 
(FWZI < 3000 km s~') whose intensity is generally 10% of the 
emission peak (GVIG99). Broad wings were observed in the 
Hell 1640 line during the 1981 bright phase. 

The optical spectra contain the weak forbidden [OIII] 4363 
line which has a relatively low critical density for collisional de- 
excitation (~10 5 cm 4 ) for 2+ . Moreover, the [OIII] 4363 line 
is observed even if [OIII] 5007, 4959, which are generally the 
strongest lines, are absent. We obtain therefore an upper limit 
on the density. The observed [OIII] 4363 line, although very 
weak, indicates a stratification of temperatures and densities of 
the emitting gas, in agreement with the profile of the physical 
conditions downstream of a shock-front. 

The Hell 4686/HeI 5876 line ratio gives information about 
the photoionization flux and the recombination conditions of 
the gas. 

In the UV, a strong NV 1240 depends on both the shock ve- 
locity and the intensity of the primary flux, while Mgll 2800, 
which is a relatively strong low-level line depends indirectly on 
the density, which speeds up the cooling rate throughout the 
nebula. Actually, Mgll can be strong if the emitting nebula 
is geometrically thick and radiation bound. The binary sepa- 
ration constrains the geometrical thickness of the radiative neb- 
ula, yielding a lower limit on the density downstream. The Hell 



3 



mo m m im ?m -r- 




iyjX' «aO Sffidu gsmo itsoo 




JD - 244000 JD - 244000 



Figure 1: Top diagram: the U and V light curves of AG Dra from 1977 adapted from Skopal et al. (2009, fig. 1, left diagram). Bottom diagram: variation of the 
strongest optical lines adapted from GVIG99, fig. 4. The red triangles indicate the JD selected for modeling (see text). 
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1640/OIII] 1660 line ratio plays a key role in the choice of the 
parameters, the density in particular. 

We have run a grid of models for each spectrum, first se- 
lecting the physical parameters which best reproduce the line 
ratios, then changing the relative abundances as a last refined 
improvement, particularly for CIV and Mgll which are the only 
representative lines of C and Mg, respectively. 

We reproduced the observed line ratios as close as possible 
considering observational errors and modeling approximation. 
A fit within a factor of 2 is generally accepted. However, at 
some epochs, the optical lines, particularly [OIII] 4363 and Hel 
4471, are underestimated by the models which explain the line 
ratios in the UV. This indicates that the emission from the ex- 
panding nebula, which is characterized by relatively low densi- 
ties, must be added. Besides stronger [OIII] and Hel lines, this 
nebula emits strong 01 and Mgll because it is extended enough 
to contain a region of recombined gas. 

After a first rough modeling of the different spectra we have 
realized that the temperature of the WD is rather stable at ~1.6x 
10 5 K on the days considered, so we iterated the modeling pro- 
cess adopting for all the models this temperature in agreement 
with Skopal et al. (2009) and Greiner et al. (1997) results 
by a fit to ROSAT quiescence data. Piro et al. (1985) de- 
rived a bb temperature of 200000 K from June 1985 EXOSAT 
observations and the Zanstra temperatures derived from the 
HeII/F(1340 A) was 109000+5400 K in quiescence (GVIG99). 

In the phases close to outbursts, both the Hell 1640/Hy6 and 
Hell 4686/Hjfi line ratios are very high. However, Hell 
1640/HyS = 14.9 cannot be reproduced by the coefficients (re- 
combination, dielectric recombination, etc.) of He and H 
adopted in the calculations (Osterbrock 1988 and references 
therein). Notice that these lines are saturated during the out- 
bursts so it is difficult to determine the broad and narrow line 
relative intensity. The best fit to the line spectra is obtained 
by adopting a relatively high shock velocity which leads to ex- 
tended regions of gas ionized to high levels. A high ionization 
parameter produces higher ratios to HB of intermediate ioniza- 
tion level lines (e.g. [OIII]/H/3 ), while a high pre-shock density 
increases the cooling rate downstream changing the stratifica- 
tion of the ions and enhancing the low ionization level lines. 

We present in Table 1 the results of modeling and we com- 
pare them with the observational data. The dates correspond- 
ing to the UV spectra are shown in column 1, followed by 
the observed UV and optical line ratios to HB . The observed 
line and continuum fluxes were reported and corrected adopt- 
ing E B _y=0.05, by GVIG99. In Table 1, next to each row con- 
taining the observed line ratios the calculated ones are shown. 
Models ml-m8 which provide the best fit to the data appear 
in the first column of Table 1 . Model m exp corresponds to the 
expanding nebula. 

Models ml-m8 and m exp are described in Table 2, in terms 
of the physical conditions in the emitting nebulae, selected as 
input parameters in the SUMA calculation process. They lead 
to the best approximation of calculated to observed line ratios. 
The first column of Table 2 shows the parameters which vary 
from epoch to epoch, while those which were found roughly 
constant are described in the bottom of the table. The shock 



velocity V s , the pre-shock density no, and the pre-shock mag- 
netic field B refer to the shock, while the color temperature of 
the hot star T, and the ionization parameter U refer to the pho- 
toionizing flux from the WD. The geometrical thickness of the 
emitting nebula D determines whether the model is radiation or 
matter bound. The dust-to-gas ratio (d/g) is crucial to explain 
the SED in the IR. The relative abundances to H of the elements 
appear in rows 6-9 of Table 2. The absolute H/3 flux calculated 
at the nebula is also given in Table 2, because it allows to cal- 
culate the absolute fluxes of all the lines. A magnetic field B () = 
10~ 3 gauss is adopted for all models. 

The calculated line ratios presented in Table 1 refer to the 
reverse shock between the stars (Sect. 2). Therefore, the shock- 
front edge of the emitting nebula downstream corresponds to 
the edge of the nebula reached by the photoionization flux from 
the WD. The velocity of the reverse shock is generally consid- 
ered approximately stable. A previous analysis (e.g. Contini & 
Formiggini 2003) concludes that the reverse shock between the 
stars is actually a standing shock, although disrupted by R-T 
and K-H instabilities. 

The values of the input parameters which appear in columns 
2-12 were selected because they explain the observational data. 
They will be considered as results of our modeling. 

In order to understand the trend of the line ratios, we show 
the profile of the electron temperature T e and electron density 
N e downstream, as well as that of the most significant ions, as 
a function of the distance from the shock-front in Fig. 2 for the 
nebulae relative to models ml, m5, and m5b (see Sect. 2.3). 
Recall that the density downstream is higher than the pre-shock 
density by a factor which depends on V s (besides the factor of 
4 determined by the adiabatic jump at the shock-front) and Bo. 
The physical conditions downstream follow the recombination 
and cooling rates, so lines from different ionization levels are 
emitted from the gas at different distances from the shock-front. 

The temperature of the gas is maximum in the immediate 
post-shock region T-1.5 10 5 ((V. s /(100 kms" 1 )) 2 , then it de- 
creases following the cooling rate determined in each slab by 
free-free, free-bound and line emission rates. At a temperature 
< 10 s K recombination produces strong line emission. How- 
ever, the temperature is maintained at < 10 4 K by the primary 
radiation from the WD and the secondary radiation from hot 
gas slabs, collisionally heated by the shock. Low ionization and 
neutral lines emitted from this region weaken in geometrically 
thin clouds. 

3.2. Broad lines 

GVIG99 present the broad component flux of Hell 4686 and 
HB corrected for extinction in their table 2. We have constrained 
the model assuming V s between 700 and 3000 km s"\ n > 10 7 
cm 3 , U is derived from the Hell h fHB h ratio, and the radius of 
the nebula within the AG Dra system (Sect. 4.2). We have 
found that V s =1000 kms" 1 , n =3 10 7 cnr 3 , and f/=100 gives 
HeII 6 /H* = 1 . The observed ratio is 1 .6. The HB"/HB h line ratio, 
calculated by this model, results in 15.2 at Earth in agreement 
with GVIG99 observations. Models m2 b , m5/,, etc., are calcu- 
lated adopting solar abundances because sputtering destroys the 
grain at such high V s (Fig. 2). 
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Table 1: Comparison of calculated with observed line ratios 1 to H/3 =1 





NV 


OI 


OIV] + 


NIV] 


CIV 


Hell 


OIII] 


Mali 


Hell 


Hy 


fOIIIl 


Hel 


Helln 




1240 


1340 


1400 


1486 


1550 


1640 


1663 


2800 


3200 


4340 


4363 


4471 


4686 


4/11/85 


1.16 


0.165 


1.06 


0.34 


1.14 


3.49 


0.34 


0.33 


0.34 


0.36 


0.001 


0.025 


0.74 


ml 


0.9 


0.18 


0.8+ 


0.6 


1.47 


4.7 


0.4 


0.32 


0.28 


0.37 


0.0014 


0.026 


0.65 


18/1/86 


2.078 


0. 


0.88 


0.37 


2.95 


14.66 


0.26 


0.2 


0.65 


0.51 


0.009 


0.049 


1.4 


ml 


2.56 


0. 


1.1 


0.23 


4.3 


7. 


0.2 


0.07 


0.4 


0.45 


0. 


0.07 


0.85 


m2,, 


0.4 


0. 


0. 


0. 


0.2 


8.9 


0. 


0. 


0.53 


0.49 


0. 


0. 


0.1 


lll exp 


0.002 


0.15 


0.02 


0.02 


0.064 


1.88 


0.23 


2.44 


0.11 


0.4 


0.04 


0.11 


0.27 


30/4/86 


2.39 


0.2 


1.26 


0.28 


2.26 


4. 


0.253 


0.0 


0.33 


0.44 


0.006 


0.045 


0.64 


m3 


2.2 


0.1 


1.8 


0.12 


1.9 


5.8 


0.26 


0.24 


0.35 


0.43 


le-4 


0.01 


0.78 


27/2/93 


0.78 


0.001 


0.35 




0.58 


2.97 


0.175 


0.184 


0.44 


0.35 


0.001 


0.036 


0.526 


m4 


0.6 


0.001 


0.32 


0.035 


0.69 


3. 


0.18 


0.17 


0.20 


0.46 


0. 


0.03 


0.42 


4/7/94 


3.9 


0.0 


3.35 


0.72 


2. 


10. 


0.63 


0.112 


0.473 


0.38 


0. 


0.093 


0.86 


m5 


3.63 


0.0 


3.7 


0.5 


3. 


6.6 


0.5 


0.12 


0.40 


0.43 


4e-4 


0.004 


0.86 


m5 fe 


0.2 


0.0 


0.17 


0.04 


16.7 


5.5 


0.2 


2.5 


0.33 


0.45 


l.e-4 


0.015 


0.75 


6/12/94 


1.24 


0.125 


1.474 


0.161 


1.55 


6.95 


0.4 


0.13 


0.35 


0.39 


0.011 


0.049 


0.64 


m6 


1.1 


0.01 


1.6 


0.35 


1.53 


6. 


0.7 


0.13 


0.3 


0.44 


0.001 


0.01 


0.78 


m6 fe 


0.06 


0.005 


0.05 


0.014 


6. 


3.7 


0.11 


4.6 


0.22 


0.46 


0. 


0.026 


0.53 


28/7/95 


1.5 


0.42 


1.61 


0.58 


2.73 


8.85 


0.5 


0.185 


0.274 


0.42 


0.017 


0.06 


0.63 


m7 


1.88 


0. 


2. 


1.1 


2.1 


7. 


0.6 


0.19 


0.42 


0.44 


0.004 


0.004 


0.89 


14/2/96 


0.52 


0.1 


0.66 


0.48 


1.14 


3.57 


0.376 


0.086 


0.12 


0.34 


0.015 


0.049 


0.53 


m8 


0.5 


0.1 


0.7 


0.46 


1.58 


4.6 


0.46 


0.054 


0.27 


0.39 


0.02 


0.03 


0.64 



1 Gonzalez-Riestra et al. (1999, Tables 1 and 2) 



Table 2: Description of the models 



model 


ml 


m2 




m3 


m4 


m5 


m5 fe 


m6 


m6/, 


m7 


ml h 


m8 




v s 


130 


100 


1000 


150 


100 


140 


700 


120 


700 


180 


700 


120 


100 


no 


5.e8 


2.4el0 


3.e7 


8.e9 


3.e9 


2.e9 


5.e8 


l.e9 


5.e8 


5.e7 


5.e8 


3.e7 


3.e6 


U 


2.5 


1 


100 


0.1 


0.5 


0.2 


5 


0.06 


5 


0.4 


5 


2 


0.1 


D 


2.6el3 


7.8el0 


4.2el3 


5.5e7 


6.ell 


8.8e8 


4.5ell 


6.e8 


4.6ell 


3.7el0 


4.6ell 


4.1el4 


l.el4 


HyS 


2.2e7 


2.8e8 


2.2e4 


1.5e5 


9.2e7 


2.8e4 


9.4e6 


2.5e4 


3.0e7 


2.6e3 


3.0e7 


3.8e6 


5.9e3 


C/H 


1.3e-5 


8.0e-6 


3.3e-4 


1.3e-5 


1.2e-5 


1.3e-5 


3.3e-4 


1.7e-5 


3.3e-4 


l.e-5 


3.3e-4 


1.3e-5 


1.3e-5 


N/H 


4.e-5 


2.1e-5 


9.1e-5 


2.1e-5 


8.1e-5 


2.1e-5 


9.1e-5 


2.1e-5 


9.1e-5 


2.1e-5 


9.1e-5 


2.1e-5 


2.1e-5 


O/H 


2.e-4 


1.6e-4 


6.6e-4 


1.6e-4 


1.6e-4 


1.3e-4 


6.6e-4 


1.6e-4 


6.6e-4 


l.e-4 


6.6e-4 


1.6e-4 


1.6e-4 


Mg/H 


2.e-5 


2.e-5 


3.e-5 


3.e-6 


2.0e-6 


1.5e-5 


3.e-5 


3.e-6 


3.e-5 


4.e-5 


3.e-5 


3.e-6 


3.e-6 


r 


8.4el3 


2.4el3 


4.2el2 


7.5el4 


2.6el3 


9.4el4 


1.8ell 


1.7el5 


1.88ell 


1.9el5 


1.3ell 


1.7el4 


2.4el4 


R 


5.8el2 


1.6el2 


6.7el2 


3.4el2 


4.8el2 


5.6el2 


4.ell 


1.6el3 


2.0ell 


2.el3 


5.3ell 


2.8el3 


2.el4 


n 


-15.8 


-17. 


-18.5 


-14. 


-16.9 


-13.5 


-21.3 


-13. 


-21.2 


-12.9 


-21.5 


-15.3 


-15. 



V s is in km s 1 , no in cm 3 , D in cm, and H/? is in erg cm 2 s 1 . 

For all the models B =10- 3 gauss, T« =160000 K, He/H = 0.1, Ne/H= 8.3 10~ 5 , Si/H=3.3 10~ 6 , S/H= 1.6 10~ 6 , A/H = 6.3 10~ 6 , 
Fe/H = 6.2 10 6 , and d/g=lQ- u by number. 
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Figure 2: Top panels: the profiles of the electron temperature, electron density and the dust temperature in the regions downstream emitting the different lines for 
models ml (left diagram), m5 (middle) and the model calculated with V s =700 km s (right diagram). The dotted red line refers to the grain radius a ?r in A. Bottom 
panels: the profiles of the most significant ion relative fractions. 



In the following active phase, starting at day 4/7/1994, we 
found that a model with V s =700 kms" 1 , no=5.10 8 cm 4 and 
U=5 reproduces the observed Uell b /H/3 b -0.8 for £>= 4.55 10 11 
cm and ~0.5 for D- 4.6 10 11 cm, as observed on 6/12/1994, and 
28/7/95, respectively. The Hell/HyS ratios in fact decrease with 
temperature at larger geometrical thickness. We do not have 
enough data to further constrain the input parameters of mod- 
els m6i, and m7^ and they are therefore identical. These mod- 
els contribute a weak broad component to the CIV line profile, 
which is not confirmed by the observations (GVIG99). Carbon 
could be included in small graphite grains which are, however, 
completely sputtered at such relatively high shock velocities. 
Therefore, depletion of carbon could be due to the presence of 
CO molecules, the second most abundant species in the ISM 
after H 2 . 

3.3. The continuum SED 

In Fig. 3 we compare the observed continuum SED with 
model calculations at each of the selected epochs. 

The analysis of the line spectra at each epoch determines the 
bb temperature of the WD and the characteristics of the nebu- 
lae downstream of the network shock-fronts. The analysis of 
the continuum SED reveals the presence of dust shells. The 
temperature of the RG is approximated by a bb flux. We use 
an effective temperature of the cool star T e //=4300 K consis- 
tent with AG Dra's classification as an early K giant (Lutz et al. 
1987) 

To constrain the continuum we have used the two data ob- 
served at 1340 and 2880 A by GVIG99. Indeed two points 
are too few to give a somehow complete information, however, 
the two wavelengths reported are strategic enough because one 
(2880 A) belongs to the bremsstrahlung from the nebulae and 
the other (1340 A) often to the bb flux from the WD. 

The nebulae are located in the circumstellar region between 
the stars and/or in the circumbinary regions. The radius r of the 
nebulae can be obtained by comparing the calculated with the 



observed absolute fluxes. They are given in the bottom of Table 
2 providing a picture of the system on a large scale (Sect. 4.2). 

In Fig. 3 we present the observations obtained in the X-ray 
and UV ranges by XMM and FUSE, respectively by Skopal et 
al. (2009, fig. 2, left panels). The data in the radio come from 
Mikolajewska (2002). We have inserted them in the spectrum 
referring to 30/4/1986, even if the epochs were not exactly the 
same, to constrain the bremsstrahlung in the radio range. Ac- 
tually, Angeloni et al. (2010) claim that radio variability also 
throughout a relatively long period does not affect the main re- 
sults of modeling. The optical low-resolution spectrum comes 
from the Loiano observatory and the flux points are determined 
by broad band optical UB VRcLc and near-IR JHKLM photom- 
etry (Skopal et al. 2009). 

4. Results 

We have modeled the line and continuum spectra emitted 
from AG Dra at different phases at epochs corresponding to 
close UV and optical observation days. The results are pre- 
sented in the following. 

4.1. Gas physical conditions 

The input parameters V s , no, and U which provide the best 
fit of the line ratios at different epochs are shown in Fig. 4 (left 
diagram). 

We found low- velocity (lv) shocks with V s between 100 and 
180 km s _I , in agreement with the FWHM of the UV high-level 
line profiles (v=100-160km s 1 ), e.g. MgVII 2510, 2629, SiVII 
2147, etc., reported by Young et al. (2006) from recent Hubble 
Space Telescope STIS observations. 

The high velocity (hv) shocks (V s = 700-1000 kms 1 ) ex- 
plain the broad wings in the Balmer lines observed close to the 
outburst epochs. They are essential in modeling the high He 
II 1640/HjS line ratios at the dates 18/1/1986 , 4/7/1994, and 
28/7/1995. The hv shocks accompany the blast wave from the 
explosion, which propagates radially, except towards the RG, 
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Figure 3: The modeling at different phases. The observations by Skopal et al (2009) are indicated by triangles: open blue at burst, filled black at quiescence, open 
black at transition. Black filled squares are the continuum data from GVIG99. Red solid line: the bb flux from the RG; solid magenta line: the bb flux from the 
WD; solid black line: the two curves represent the bremsstrahlung from the nebulae described in Table 2 and reprocessed radiation by dust; long-dashed black line: 
the bremsstrahlung from the high velocity models; short-dashed black line: the expanding nebula; dotted line: bremsstrahlung from a model calculated by V s =280 
km s~' . 
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where collision with the dense RG wind prevents its expansion 
and creates a complex network of shock-fronts. 

The SEDs of the continuum emitted from the nebulae down- 
stream of the hv shock-fronts are very different from those 
reported by Skopal et al. in quiescent and transition phases 
(Fig. 3). In particular they produce negligibly low UV-optical 
fluxes. On the other hand the continuum SEDs emitted from 
the nebulae downstream of the lv shock-fronts between the 
stars can reproduce the trend of the observations in the soft X 
ray range for V s > 150kms _1 . For V s > 180 km s~' the X-rays 
can even overcome the bb flux from the WD, mimicking a 
higher WD temperature. For example, the SED calculated with 
V s =280 kms^'in the transition phase on 27/2/1993 is shown 
in Fig. 3. 

The pre-shock densities are relatively high, leading to the 
low velocities accompanying the reverse shock between the 
stars in AG Dra (from conservation of mass at the shock-front: 
n0v0=nlvl, where refers to upstream and 1 to downstream 
conditions). In R Aqr the velocities of the reverse shock 
are ~100 km s -1 and the densities are ~10 5 cirT 3 (Contini & 
Formiggini 2003). The velocities of the reverse shock inferred 
from the IR line ratios in D-type SS (Angeloni et al. 2007a) 
are ~400-500 kms -1 . We wonder whether dust and/or higher 
densities could be responsible for the deceleration of the 
reverse shock. 

The input parameter U varies between a maximum of 100 in 
the high velocity nebulae to a minimum of 0.06 for model m6. 
A high U indicates that the nebula is close to the hot source 
and/or that there is no matter obstructing the photoionizing flux 
(cfr. for Z And - Contini et al. 2010, in preparation). 

4.2. Radius of the nebulae 

The comparison of the continuum fluxes calculated at the 
nebula with those observed at Earth determines the radius r of 
the nebula by the factor rj (r 2 = XO'* 1 d 2 ff, where ff is the filling 
factor, d the distance to Earth of AG Dra). Adopting d=2.5 Kpc 
(GVIG99) and |=lwe obtain the radius upper limits which 
appear in row 10 of Table 2. Considering that the binary sepa- 
ration is 2.3 10 13 cm (Leibowitz et al. 1985, Tomov & Tomova 
2002) ff is likely to be -0.1- 0.01. 

We have calculated the distance of the inner edge of the nebu- 
lae from the WD by combining F (the number of photons cm -2 
s _1 reaching the nebula) with the ionization parameter U calcu- 
lated phenomenologically by modeling the line spectra, 
F (r WD /R) 2 = Unc, where r WD is the radius of the WD (0.08 R 
GVIG99), n the density and c the speed of light. The resulting 
R are given in row 1 1 of Table 2. 

Model (m exp ) should be added to the reverse shock models at 
the dates 30/4/1986, 4/7/1994/, and 28/7/1996, increasing the 
flux of the low-ionization and neutral lines. This model repre- 
sents the nebula downstream of the shock propagating out of 
the system (Contini et al. 2009c). The shock-front edge of the 
nebula is opposite the edge reached by the flux from the WD. 
no is lower and the nebula is extended. On the other hand, the 



nebulae downstream of the hv shocks (700-1000 km s~') have 
a very small radius because they represent matter ejected at the 
outburst which had not yet the time to propagate in the sur- 
roundings of the WD. 

4.3. Integrated UV and soft X-ray fluxes from the nebulae 

The input parameters (Fig. 4 left) explain the trend of the 
integrated fluxes (Fig. 4 right) in the optical F op , (8267-3100 
A), UV F uv (3 100-620 A), and in the soft X-ray range F x (0. 1 
-1.2 keV). The maximum gas densities downstream, i.e. the 
densities which are used to calculate the fluxes for each model 
are shown in the right diagram, as well as the E.M. (the emis- 
sion measure ZDn n e ). It can be noticed that F op , and F^y 
follow the trend of U, while F^ increases with V s , because the 
maximum temperature downstream is oc V 2 . 

A rough anti-correlation (Fig. 5) is evident between Fx and 
Fuv calculated at the nebula by models ml - m8, m2/,, and m5i, 
. The fluxes calculated in models m6b and m7 J >, coincide with 
m5b, therefore, they are removed from Fig. 5 for the sake of 
clarity. Recall that the higher V s the larger is the high tempera- 
ture zone within the nebula because recombination coefficients 
are low. The cooling rate increases rapidly when the tempera- 
ture of the nebula decreases below 10 5 K. At this stage a strong 
compression followed by a rapid cooling rate (oc n 2 ) reduces 
the geometrical thickness of the UV emitting region (Fig. 2, 
diagram on the right) and the UV flux is relatively low. So the 
SED of the continuum from hv and lv shocks are different. In 
the soft X-ray range the flux from lv shocked nebulae dominates 
over the WD bb only for V s > 150 km s" 1 . 

The comparison of the fluxes calculated at the nebula with 
the continuum SED observed at Earth indicates that hv shocks 
have a radius smaller than that of lv ones (Table 2, last row), 
so the fluxes observed at Earth of hv nebulae would be lower 
by several orders of magnitude than those corresponding to lv 
nebulae and WD bb radiation. However, both the latter fluxes 
are absorbed by the envelop of debris accompanying the explo- 
sion. The drop in the X-ray counts at the outbursts can thus be 
explained by the different origin of the X-ray. 

Summarizing, at outbursts, the X-rays are emitted from the 
shock accompanying the blast wave (which absorbs the WD bb 
flux). The observed flux is relatively low because it is diluted by 
the small rj factor. In the quiescent and transition phases the soft 
X-ray data (observed at similar phases but at different epochs by 
Skopal et al) are reproduced by the WD bb with a contribution 
from the bremsstrahlung which dominates for shock velocities 
> 150 kms~'. The optical-UV flux is bremsstrahlung emitted 
from downstream nebulae photoionized by the WD bb flux. 

4.4. Relative abundances 

The abundances of the heavy elements relative to H which 
explain the observed line ratios are given in Table 2 (rows 6-9). 
We could model only the elements corresponding to observed 
lines. 

The results show a strong under-abundance of the heavy ele- 
ments relative to the solar ones (C/H=3.3 10~ 4 , N/H=9.1 10~ 5 , 
0/H=6.6 10" 4 , Mg/H=2.6 10" 5 , Nussbaumeret al. 1988). 
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Figure 4: The input parameters (left diagram), the post-shock density, the fluxes in the optical, UV, and X-ray domain, and E.M., calculated by the models (right 
diagram) 
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Figure 5: F^ versus F UB calculated at the nebula by models which explain the 
line and continuum spectra. Filled triangles: low-velocity models; open circles: 
high- velocity models. 



Nussbaumer et al. (1988) compared the O/N and C/N ratios 
for different objects, e.g. symbiotic stars, planetary nebulae, 
novae, etc. Our results show that AG Dra has C/N compara- 
ble with other symbiotic stars, but higher O/N. Moreover O/N 
ratios are even outside the ranges covered by the other objects. 

We find C/O varying from 0.05 to 0.13, while the solar abun- 
dance ratio is 0.5. The N/H ratio is underpredicted by a factor 
of 4 relative to the solar, while C/H is lower than solar by a fac- 
tor > 10. Also O/H is lower than the solar by a factor of ~2.7 
in average, while Mg/H varies between solar and 0. 1 solar de- 
pending on the phase. Mg can be trapped into silicate grains, so 
the different abundance ratios indicate an inhomogeneous dis- 
tribution of dust. We cannot comment about the Si/H relative 
abundance, because no Si lines are given in the spectra. SilV] 
line is blended with OIV] in the 1400 A multiplet. The calcu- 
lated spectra show that the intensity of the SilV] line is always 
lower than that of OIV] or similar in the case of model m4. For 
m4 N/H is rather solar, while C and O are depleted strengthen- 
ing the hypothesis that C and O are trapped into CO molecules. 

The origin of the high depletion of C, N, and O relatively to H 
and He is most probably due to inclusion in dust grains and/or 
molecules (Angeloni et al. 2007a). Dust grains are coupled 
to the gas by the magnetic field downstream of shock-fronts. 
They are heated by radiation from the WD and collisionally 
by the gas. Starting by an average grain radius a gr =0.2 /mi, 
Fig. 2 (top panels) shows that the grains reach a temperature 
< 2000 K and survive sputtering and evaporation (which are 
calculated consistently with the temperature of the gas) in the 
transition and quiescence phases, but are completely destroyed 
downstream of high velocity shocks. 

Fig 3 shows that dust reprocessed radiation is hidden 
throughout the SED by the bb flux from the RG. Adopting dig 
ratios (Table 2) similar to those of the ISM (4.1 10~ 4 by mass 
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which corresponds to 10 by number) dust emission does not 
contribute to the continuum sensibly. At the dates 6/12/1994 
and 28/7/1995 a d/g ~10~ 15 has been adopted. However, these 
results are based on the ISO observations at date 29/July/1996. 
Actually, dust emission can vary with phase, as was found for 
e.g. RS Ophiuci by Rushton et al. (2010). 

Oxygen is depleted because it is trapped in unsputtered sil- 
icate grains and in CO molecules, while carbon is much more 
depleted (by a factor of 3) because it is trapped into graphite 
grains, PAHs, and diatomic molecules like CO, CN, CS, etc. 
(Contini & Shaviv 1982). Formation of dust species down- 
stream of shocks is characteristic of the RG. Actually, shock- 
excited far-IR emission of CO was detected by Reach & Rho 
(1998) with the ISO Long-Wavelength Spectrometer from the 
SN remnant 3C391. 

As a result, the depletion of C, N, O and the ~1 year fluctu- 
ation period observed during the high outbursts in the U band 
and in the V band strengthens the hypothesis that this short pe- 
riod depends on the sudden disturbance of the dynamical field 
between the stars by the collision of the dusty shells ejected 
from the RG with the wind from the WD. 

5. Discussion and concluding remarks 

In previous sections we have demonstrated that the 
bremsstrahlung emitted downstream of shocked nebulae in AG 
Dra SS contributes to the soft X-rays at different phases. 

GVIG99 report that after the quiescence period from 1969 to 
1980 and the active phase from November 1980 and beginning 
1983, a new minor outburst took place in March 1985 followed 
by another one after 340 days in January 1986. After 1986 a 
8 years period of quiescence followed until the major active 
phase started in June 1994 with at least five peaks in June 1994, 
June 1995, July 1996, and June 1997, and one in summer 1998. 
They are separated by nearly the same time interval of about 
360 days. Another small light maximum is present 330 days 
before the first maximum (Fig. 1, top). 

We have reproduced the UV-optical spectra at different 
epochs by model calculations (Tables 1 and 2): Nov 4 1985 
(ml), Jan 18 1986 (m2), Apr 30 1986 (m3), Feb 27 1993 (m4), 
Jul 4 1994 (m5), Dec 6 1994 (m6), Jul 28 1995 (m7), and Feb 
14 1986 (m8). The bremsstrahlung calculated by model ml 
refers to a transient phase, m2 to a maximum, m3 to the end of 
the maximum, near quiescence, m4 to quiescence, m5 is close 
to maximum, m6 corresponds to a transient phase in terms of 
a small depression between two maxima, m7 is close to maxi- 
mum, and m8 represents a quiescence phase. 

The classification in quiescence, maximum, and transient 
phases by Skopal et al, provides only a schematic picture of 
the physical conditions. We believe that the conditions change 
from burst to burst, from quiescence to quiescence, etc., be- 
cause collision of the fast wind from the WD with the RG wind 
between the stars and with the medium inhomogeneities on the 
side opposite the RG, at outbursts, leads to a network of shock- 
fronts which varies from time to time. 

Accepting that Skopal et al. (2009, fig. 2) continuum data 
can be applied at different epochs, Fig. 3 shows that the soft 



X-rays would correspond to the bb flux from the WD in tran- 
sient and quiescence epochs for shocks with V s < 150 kms~'. 
During the 18/1/1986 burst the X-rays originated from the high 
velocity shock and could be observed also beyond 1 keV. After 
the burst, on 30/4/1986, the soft X-rays corresponding to the bb 
flux are comparable with the bremsstrahlung downstream of the 
reverse shock with V s =150 kms~'. 

During the 4/7/1994 and 28/7/1995 bursts the bremsstrahlung 
downstream of the reverse shock corresponds to the soft X-ray 
range, while downstream of the blast wave shock it covers also 
the X-ray range (> 1 keV). There are no contemporary data 
which could confirm our prediction in the X-ray range. 

The fluctuations in the U band during the active phases be- 
tween 1995 and 1999 show a ~360 days periodicity (Fig. 1, 
top diagram) which is similar to the red giant pulsation period 
(~1 year). Pulsations provoke dusty shells ejection. Collisions 
of the shells with the WD wind give origin to the sequence of 
peaks during the active phases. This is connected to the gas and 
dust composition in shocked nebulae. Calculation results show 
a depletion of C by a factor > 10, a depletion of O by a factor 
of 2.7 and of N by a factor of ~4, while Mg is more variable. 
The depleted elements are trapped into dust grains (silicates and 
graphite) and/or into diatomic molecules (CO, CN, CS, etc.). 

The prerogative of observing the effect of the RG pulses 
throughout the light curve in U belongs mainly to systems 
which are seen nearly face-on, like AG Dra. Most informa- 
tion is lost in edge-on systems due to eclipses not only by the 
two main stars, but also by the different nebulae and shells (e.g. 
CH Cyg, Contini et al. 2009c) within the system. 

Fig. 3 diagrams show that the continuum in the optical-UV 
range is emitted by the shocked nebulae between the stars. The 
low-frequency tail of the WD bb flux is seen at v > 10 15 Hz. 

It was concluded by Skopal et al. (2009) that the super soft 
X-ray bump represents directly the WD bb flux. The peak of 
the WD bb between ~3 10 15 and 3 10 16 Hz is not reproduced 
by shock dominated models. Actually, it was neither confirmed 
by the observations due to the strong absorption in this range. 
Moreover, the H column density N# ~3 10 20 cirT 2 adopted by 
Skopal et al. (2009) to correct the data to fit the Planck func- 
tion, can sensibly vary with time in a complex system such as 
a symbiotic one, if column densities from the emitting nebulae 
are accounted for. 

The anti-correlation between the X-ray and UV fluxes, par- 
ticularly at outbursts, was explained by the obscuration of the 
WD by the envelop of matter ejected at outburst which absorbs 
the soft X-ray emission (Skopal et al. 2009). We suggest that 
the absorbing envelop corresponds to the fragmented nebula 
downstream of the high velocity shock accompanying the blast 
wave, which was successfully invoked to explain some SS fea- 
tures, such as the broad Lya line in CH Cyg (Contini et al. 
2009a). The fast debris emit the broad line spectra (Sect. 3.2) 
and the X-ray bremsstrahlung observed at outbursts. 

Finally, the results presented for AG Dra are in agreement 
with Smith et al. (1996) who claim by simple arguments based 
on the luminosity function of solar and low-metallicity K giants 
that yellow symbiotic stars arise among low-metallicity stars, 
because low-metallicity K giants have a larger mass loss rate, 
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thus triggering symbiotic-like activity in a binary system. Gen- 
erally yellow SS show IR emission from dust (e.g. HD 330036, 
Angeloni et al. 2007b) which cannot be observed in AG Dra 
because it is hidden the bb RG flux throughout the SED. How- 
ever, we have found that dust is present between the stars with 
a non-negligible dust-to-gas ratio (< 4 10~ 4 by mass). 
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